matlab循环修改变量名,修改MATLAB语句，实现两个变量的循环

clear all;clc;

global hx0 hy0 hz0         % 外加磁场参数

hz0=0; hy0=0;

i=0;

mz3=0;

mz4=0;

for hx0=0:0.01:5;

y0=[0.11;0.01;pi-0.11;pi+0.01];

t0=0:0.1:500;

[t,y]=ode45('rhssys',t0,y0);

mx1=sin(y(:,1)).*cos(y(:,2));

my1=sin(y(:,1)).*sin(y(:,2));

mz1=cos(y(:,1));

mx2=sin(y(:,3)).*cos(y(:,4));

my2=sin(y(:,3)).*sin(y(:,4));

mz2=cos(y(:,3));

i=i+1;

mz3(i)=cos(y(5000,1));

mz4(i)=cos(y(5000,3));

end

hx0=0:0.01:5;

subplot(2,1,1);

plot(hx0,mz3,'LineWidth',2);

xlabel('h_x0','fontsize',20);ylabel('m_z1','fontsize',20);

title('eta=0.33','fontsize',20);

set(gca,'FontSize',20);

subplot(2,1,2);

plot(hx0,mz4,'LineWidth',2);

xlabel('h_x0','fontsize',20);ylabel('m_z2','fontsize',20);

set(gca,'FontSize',20);

此循环程序现在实现了hx0和t0的双循环吗?即每一个hx0，此程序算了t0从0到500所对应的mx,my,mz吗？如果没有，如何实现hx0和t0的双循环?最后怎样改变语句，可以读出所有t0对应的y0值，而不是只是最后一组y0.